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Abstract 

In light of the recent LHC Higgs search data, we investigate the pair production of a SM-like 
Higgs boson around 125 GeV in the MSSM and NMSSM. We first scan the parameter space of each 
model by considering various experimental constraints, and then calculate the Higgs pair production 
rate in the allowed parameter space. We find that in most cases the dominant contribution to 
the Higgs pair production comes from the gluon fusion process and the production rate can be 
greatly enhanced, maximally 10 times larger than the SM prediction (even for a TeV-scale stop 
the production rate can still be enhanced by a factor of 1.3). We also calculate the x 2 value with 
the current Higgs data and find that in the most favored parameter region the production rate 
is enhanced by a factor of 1.45 in the MSSM, while in the NMSSM the production rate can be 
enhanced or suppressed {osusy I '&sm varies from 0.7 to 2.4). 

PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 
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I. INTRODUCTION 



Based on the combined data collected at the center-of-mass energies of 7 TeV and 8 TeV, 
the experimental programme to probe the mechanism of electroweak symmetry breaking 
at the LHC has recently witnessed the discovery of a new particle around 125 GeV 
The properties of this particle, according to the updated analyses of the ATLAS and CMS 
collaborations at the end of 2012 [2], roughly agree with the the Standard Model (SM) 
prediction and thus it should play a role in both the symmetry breaking and the mass 
generation. However, the issue of whether this particle is the SM Higgs boson is still open, 
and indeed there are some motivations, such as the gauge hierarchy problem and the excess in 
the di-photon channel over the SM prediction jl[ Q , to consider new physics interpretation 
of this boson. Studies in this direction have been performed intensively in low energy 
supersymmetry (SUSY) and it was found that some SUSY models can naturally provide a 
125 GeV Higgs boson 0-5], and fit the data better than the SM jf| (similar studies have also 
been performed in some non-SUSY models like the little Higgs models and two-Higgs-doublet 
or Higgs-triplet models Q). 

After the discovery of the Higgs boson, the next important task for the LHC is to test the 
property of this Higgs boson by measuring all the possible production and decay channels 
with high luminosity. Among the production channels, the Higgs pair production is a rare 
orocess at the LHC. Since it can play an important role for testing the Higgs self-couplings 
8|, |9] (the determination of the Higgs self-couplings is of great importance since it is in- 
dispensable to reconstruct the Higgs potential), it will be measured at the LHC with high 
luminosity. 

In the SM the Higgs pair production at the LHC proceeds by the parton process gg — > hh 
through the heavy quark induced box diagrams and also through the pro duction of an off- 
shell ffiggs which subsequently sphts into two ou-shell Higgs bosous Q,Q. The production 
rate is rather low for y/s = 14TeV, about 20 fb at leading order |10J] and reaching roughly 35 
fb after including the next-to-leading order QCD correction Q. The capability of the LHC 
to detect this production process was investigated in 12l-ll5l|. These analyses showed that for 
a 125 GeV Higgs boson the most efficient channel is gg — > hh — >■ 6677 with 6 signal events 
over 14 background events expected for 600 fb -1 integrated luminosity after considering 



some elaborate cuts [12| (the detection through other channels like hh — > bbW + W and 
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hh — > bbr + r has also been studied recently {13I In principle, the capability can be 



16| is applied for the 



further improved if the recently developed jet substructure technique 
Higgs tagging. 

The Higgs pair production at the LHC may also be a sensitive probe for new physics. In 
supersymmetric models such as the Minimal Supersymmetric Standard Model (MSSM) [171 ]. 
the pair production of the SM-like Higgs boson receives additional contributions from the 
loops of the third generation squarks and also from the parton process bb — > Hi — » hh with 



Hi denoting a CP-even non-standard Higgs boson [18|, [19| . It was found that in some cases 
(e.g., a light stop with a large trilinear soft breaking parameter A t and/or a large tan/3 
together with moderately light Hi), these new contributions may be far dominant over the 
SM contribution, and as a result, the rate of the pair production may be enhanced by several 
orders 
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191 ] . Note that since the experimental constraints (direct or indirect) on the SUSY 
parameter space have been becoming more and more stringent, the previous MSSM results 
should be updated by considering the latest constraints. This is one aim of this work. To 
be specific, we will consider the following new constraints: 

• The currently measured Higgs boson mass rrih = 125 GeV [2]. In SUSY this mass is 
sensitive to radiative correction and thus the third generation squark sector has been 
tightly limited. 

• The LHC search for the third generation squarks [20]. So far although the relevant 
bounds are rather weak and usually hypothesis-dependent, it becomes more and more 
clear that a stop lighter than about 200 GeV is strongly disfavored. 



The observation of B s — > by the LHCb 2l|. In the MSSM it is well known that 

the branching ratio of EL — > is proportional to tan 6 /3 /mjj for a large tan/3 and 

a moderately light H [22] . Since the experimental value of B s — > coincides well 

with the SM prediction, tan (3 as a function of has been upper bounded. 

The LHC search for a non-standard Higgs boson H through the process pp — > H — > 
t + t~ [23). Such the search relies on the enhanced Hbb coupling and the nought signal 
seen by the LHC experiments implies that a broad region in the tan /3 — m# plane has 
been ruled out. 
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The global fit of the SUSY predictions on various Higgs signals to the Higgs data re- 



ported by the ATLAS and CMS collaborations |24j, the dark matter relic density 25 1 



261 ] can also limit SUSY pa- 



as well as the XENON2012 dark matter search results 
rameters in a complex way. 

Another motivation of this work comes from the fact that the Next-to-Minimal Supersym- 



metric Standard Model (NMSSM) [27j] is found to be more favored by the Higgs data and 
the fine-tuning argument [6]. So far the studies on the Higgs pair production in the NMSSM 
are still absent. So it is necessary to extend the study to the NMSSM. 

This paper is organized as follows. In Sec. II we briefly introduce the features of the Higgs 
sector in the MSSM and NMSSM. Then in Sec. Ill we present our results for the Higgs pair 
production in both models. Some intuitive understandings on the results are also presented. 
Finally, we summarize our conclusions in Sec. IV. 



II. HIGGS SECTOR IN MSSM AND NMSSM 

As the most economical realization of SUSY in particle physics, the MSSM [17] has 
been intensively studied. However, since this model suffers from some problems such as the 
unnaturalness of fi parameter, it is well motivated to go beyond this minimal framework. 
Among the extensions of the MSSM, the NMSSM as the simplest extension by singlet 
field [271 has been paid much attention. The differences between the two models come from 
their superpotentials and soft-breaking terms, which are given by 



W MSSM = Y U Q ■ H U U - Y d Q ■ H d D - Y e L ■ H d E + ^H u ■ H d , (1) 
^nmssm = Y U Q ■ H U U - Y d Q ■ H d D - Y e L ■ H d E + \H U ■ H d S + ±kS*, (2) 
KT M = ™l\Hu\ 2 + m d \H d \ 2 + (BfiH u ■ H d + h.c), (3) 
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2 , ~ 2 

+ m d 


H d \ 


\H U 


2 + m 2 d 


H d \ 



K NMSSM = ^2 {Hul 2 + ^2 lHdl 2 + ^2 + . ^ + ^8 + ^ (4) 

Here Hi {i = u, d) and S denote gauge doublet and singlet Higgs superfields respectively, 
Q, U, I), L and E represent matter superfields with Yi (i = u, d, e) being their Yukawa 
coupling coefficients, fhi (i = u,d,S), B, A\, and A K are all soft-breaking parameters and 
the dimensionless parameters A and k reflect coupling strengthes of Higgs self interactions. 
Note the /i-term in the MSSM is replaced by Higgs self interactions in the NMSSM, so when 
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the singlet field S develops a vacuum expectation value s, an effective \x is generated by 
Mf = Xs - 

Like the general treatment of the multiple-Higgs theory, one can write the Higgs fields in 
the NMSSM as 

<j>d+i<Pd 



V u + 



V2 



H, 



v d + 



V2 



s 



h: 



(5) 



and diagonalize their mass matrices to get Higgs mass eigenstates: 
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(6) 



Here Hi, H 2 , H 3 with convention < m^ 2 < m^ 3 and Ai, A 2 with convention < 
rriAi denote the physical CP-even and CP-odd Higgs bosons respectively, G° and G + are 
Goldstone bosons eaten by Z and W bosons respectively, and H + is the physical charged 
Higgs boson. The Higgs sector in the MSSM can be treated in a similar way except that it 
predicts only two physical CP-even states and one physical CP-odd state, and consequently, 
the rotation matrices Uh and XJ a are reduced to 2 x 2 matrices. 

One distinct feature of the MSSM is that Hi usually acts as the SM-like Higgs boson 
(denoted by h hereafter) and its mass is upper bounded by mz at tree level. Obviously, to 
coincide with the LHC discovery of a 125 GeV boson, large radiative correction to nth is 
needed, which in turn usually requires the trilinear soft breaking parameter A t to be large. 
For example, in the case of large rriA and moderate tan 0, is given by 4| 

i 

,2 



i\ — M z cos 2 2/3 + 3m * 



4ir 2 v 2 



In -4 + -4 1 



mi 



12m? 



(7) 

where the first term is the tree-level mass and the last two terms are the dominant corrections 
from the top-stop sector, m 4 ~ = y/m^ nv[ 2 {m^ denotes stop mass with convention < mj 2 ) 
represents the average stop mass scale and X t = A t — \i cot 0. One can easily check that for 
a 500 GeV and 1 TeV stop, \A t \ should be respectively about 1.8 TeV and 3.5 TeV to give 
m h ~ 125 GeV. 

In the NMSSM, mh exhibits at least two new features One is that it gets additional 



contribution at tree level so that m\ tree 



AV) cos 2 2(3 + AV, and for A ~ 0.7 and 



tan f3 ~ 1, mh can reach 125 GeV even without the radiative correction. The other feature is 




FIG. 1: Feynman diagrams for the pair production of the SM-like Higgs boson via gluon fusion in 
the MSSM and NMSSM with Hi denoting a CP-even Higgs (1 = 1,2 for the MSSM and I = 1,2,3 
for the NMSSM) and qi j = 1, 2) for a squark. The diagrams with initial gluons or final Higgs 
bosons interchanged are not shown here. For the quarks and squarks we only consider the third 
generation due to their large Yukawa couplings. 

that the mixing between the doublet and singlet Higgs fields can significantly alter the mass. 
To be more explicit, if the state Hi is h, the mixing is to pull down the mass, while if H2 
acts as h, the mixing will push up the mass. Another remarkable character of the NMSSM 
is that in the limit of very small A and k (but keep \i fixed), the singlet field decouples from 
the theory so that the phenomenology of the NMSSM reduces to the MSSM. So in order to 
get a Higgs sector significantly different from the MSSM, one should consider a large A. 

Throughout this work, we require 0.50 < A < 0.7 in our discussion of the NMSSM and 
we consider two scenarios: 
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FIG. 2: Feynman diagrams for the parton process bb — > hh in the MSSM and NMSSM. 

• NMSSM1 scenario: Hi acts as the SM-like Higgs boson. For this scenario, the addi- 
tional tree-level contribution to rrih is canceled by the mixing effect, and if the mixing 
effect is dominant, the parameters in the stop sector will be tightly limited in order to 
give m h ~ 125 GeV. 

• NMSSM2 scenario: H 2 acts as the SM-like Higgs boson. In this scenario, both the 
additional tree-level contribution and the mixing effect can push up the mass. So for 
appropriate values of A and tan/3, m h can easily reach 125 GeV even without the 
radiative correction. 



III. CALCULATIONS AND NUMERICAL RESULTS 

In SUSY the pair production of the SM-like Higgs boson proceeds through the gluon 
fusion shown in FigfJJand the bb annihilation shown in Figj2j These diagrams indicate that 
the genuine SUSY contribution to the amplitude is of the same perturbation order as the 
SM contribution. So the SUSY prediction on the production rate may significantly deviate 
from the SM result. To ensure the correctness of our calculation, we checked that we can 



reproduce the SM results presented in [10|| and the MSSM results in 18]. Since the analytic 
expressions are quite lengthy, we do not present here their explicit forms. 

In our numerical calculation we take m t = 173 GeV, = 4.2 GeV, mz = 91.0 GeV, 



my/ = 80.0 GeV and a = 1/128 [28], and use CT10 [29j to generate the parton distribution 
functions with the renormalization scale /ir and the factorization scale \xf chosen to be 2rrih- 
The collision energy of the LHC is fixed to be 14 TeV. Then we find that for = 125 GeV, 
the production rate in the SM is 18.7 fb for gg — > hh and 0.02 fb for bb — > hh (the rates 
change very little when m h varies from 123 GeV to 127 GeV). 
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FIG. 3: The scatter plots of the surviving samples, showing asiJSY /&SM versus the SM-like Higgs 
boson mass. The plus '+' (blue) denote the results with only the gluon fusion contribution, while 
the circles V (pink) are for the total results. 



For each SUSY model we use the package NMSSMTools-3.2.0 30[] to scan over the pa- 
rameter space and then select the samples which give a SM-like Higgs boson in the range of 
125 ± 2 GeV and also satisfy various experimental constraints, including those listed in Sec- 
tion I. The strategy of our scan is same as in [(J except for three updates. First, since the rar e 



decay B s — > has been recently observed with Br(B s — > = 3.2l 12 x 10~ 9 21], 

we use a double-sided limit 0.8 x 10~ 9 < Br(B s — > /i + /U~) < 6.2 x 1CT 9 . Second, for the 
LHC search of the non-standard Higgs boson, we use the latest experimental data 23J]. The 
third one is that we require stops heavier than 200 GeV {20]. After the scan, we calculate 
the Higgs pair production rate in the allowed parameter space. We will demonstrate the 
ratio (Tsusy/o'sm for each surviving sample. Of course, such a ratio is less sensitive to higher 
order QCD corrections. 

In Fig. [3] we show the normalized production rate as a function of the Higgs boson mass 
for the surviving samples in the MSSM and NMSSM (for the NMSSM we show the results 
for the NMSSM1 and NMSSM2 scenarios defined in Sec. II). This figure shows two common 
features for the three scenarios. One is that the production rate can deviate significantly 
from the SM prediction: in most cases the deviation exceeds 30% and in some specail cases 
the production rate can be enhanced by one order. The other feature is that for most 
cases the dominant contribution to the pair production comes from the gluon fusion, which 
is reflected by the approximate overlap of V (pink) with '+' (blue). Fig. |3] also exhibits 
some difference between different scenarios. For example, in the MSSM the bb annihilation 
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contribution can be dominant for some surviving samples, which, however, never occurs in 
the NMSSM. Another difference is that the NMSSM1 tends to predict a larger production 
rate than other scenarios. 

Now we explain some features of the results in Fig. [3J First, we investigate the cases of 
the MSSM where the bb annihilation plays the dominant role in the production. We find 
that they are characterized by a moderately large tan /3 (tan (3 ~ 10 so that the Hbb coupling 
is enhanced), a moderately light H (300 GeV < ran ^ 400 GeV) and a relatively large Hhh 
coupling. While for the NMSSM scenarios, since we are considering large A case, only a 
relatively small tan f3 is allowed so that the Hibb coupling is never enhanced sufficiently js|. 
We also scrutinize the characters of the gluon fusion contribution in the MSSM. As the first 
step, we compare the sbottom loop contribution with the stop loop. We find that for the 
surviving samples the former is usually much smaller than the latter. Next we divide the 
amplitude of Fig. [I] into five parts with M 1? M 2 , M 3 , M4 and M 5 denoting the contributions 
from diagrams (l) + (2), (3) + (4), (5), (6)+(7) and (8) + (9) + (10), respectively. For each of 
the amplitude, it is UV finite so we can learn its relative size directly. We find that the 
magnitudes of M 2 and M3 are much larger than the others. This can be understood as 
follows: among the diagrams in Fig. [H only (3), (4) and (5) involve the chiral flipping of the 
internal stop, so in the limit m^mj 3> Im^ the main parts of M 2 and M3 can be written 
as 

M ~ a 2 Y 2 ( Cl sin 2 26^ + c 2 ^) (8) 

m~ m- 

tl t-2 

where Y t is the top quark Yukawa coupling, 9 t and A t are respectively the chiral mixing 
angle and the trilinear soft breaking parameter in the stop sector, and c\ and c 2 are 0(1) co- 
efficients with opposite signs. Since a large A t is strongly favored to predict rrih ~ 125 GeV 
in the MSSM 3[ and the other contributions are usually proportional to m 2 /m 2 or m^/m 2 , 
one can easily conclude that M 2 and M 3 should be most important among the five ampli- 
tudes. In fact, we checked that without the strong cancelation between M 2 and M3, the 
production rate can easily exceed 100 fb for most surviving samples. 

As a proof for the validity of Eq. (jSJ) , in Fig. H] we show A t /mi versus mt , where the 
samples are classified according to the value of R = &susY(gg — > hh)/asM{gg hh). The 
left panel indicates that in the MSSM the region characterized by a light m^ and a large 
|A/ m *il usually predicts a large R. This can be understood as follows. In the MSSM with 
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FIG. 4: Same as Fig[3l but showing Af/m? versus mx . The samples are classified according to 
the value of R = crsuSY(gg —> hh)/asM(gg — > hh) with a denoting the hadronic cross section via 
gg -)• hh. 

a light ti, the other stop (£2) must be sufficiently heavy in order to predict ~ 125 GeV 
[3I. Then, after expressing sin 2 26 t in terms of A t and stop masses, one can find that the first 
term in Eq.® scales like (A t /mtJ A (m^m^/m~^), and therefore its value grows rapidly with 
the increase of |A t /mjJ and is unlikely to be canceled out by the second term in Eq.fjHJ). In 
fact, the upper left region of the panel reflects such a behavior. This panel also indicates 
that even for t± and ti at TeV scale, the production rate in the MSSM may still deviate from 
its SM prediction by more than 30%. This is obvious since \A t \ in Eq.® is usually larger 
than stop masses [3J. Finally, we note that for mi > 1 TeV, there exist some cases where 
the deviation is small even for A t /mi ~ 3. We checked that these cases actually correspond 
to a small mass splitting between t\ and £2- I n such a situation, the first term in Eq.flBJ) is 
proportional to A%/m~ (since t — 7r/4), and its contribution to the rate is severely canceled 
by the second term. 

Eq.Q may also be used to explain the results of the NMSSM1 scenario. In this scenario 
we checked that the mixing effect on often exceeds the additional tree level contribution 
(as discussed in Sec. II), and consequently the soft breaking parameters in the stop sector 
are more tightly limited than the other two scenarios. For example, given the same values 
of m% x and m~ t2 for the three scenarios, the NMSSM1 scenario usually prefers a larger \A t \. 
Consequently, this scenario tends to predict the largest production rate according to Eq.(|HJ). 
As for the R value in the NMSSM2 scenario, the situation is quite complex because a large A 
alone can push the value of up to about 125 GeV and thus the soft breaking parameters 
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FIG. 5: Same as Fig. El but showing cjsusy / &sm versus \ 2 ■ Here only the samples satisfying 125 
GeV <m h < 126 GeV are plotted. 

in the stop sector are not so constrained by the Higgs mass [3j. But, anyway, this scenario 
still has the features that R is maximized for a large A t and a light t\ and that R can deviate 
sizably from unity for TeV-scale stops. 

Finally, we focus on the samples which predict a SM-like Higgs boson in the best fitted 
mass region, 125 GeV < m h < 126 GeV J24L For these samples, we calculate the x 2 value 
with the LHC Higgs data (for details, see 6|,|24j) and show its correlation with the normalized 
rate (Jsusy/o'sm in Fig- El This figure indicates that in the MSSM and NMSSM2 scenarios, 
there exist a lot of samples with x 2 much smaller than its SM value (x|m = 16.5), which 
implies that the MSSM and NMSSM2 scenarios may be favored by the current data |6j. In 
contrast, the NMSSM1 scenario can only slightly improve the fit. From this figure we also 
see that in the favored parameter space with a small x 2 the production rate can sizably 
deviate from the SM prediction ( in the parameter space with a large x 2 the production rate 
can be several times larger than the SM value). For example, in the low x 2 region of the 
MSSM, the normalized rate is approximately 1.45, while in the NMSSM2 scenario the rate 
varies from 0.7 to 2.4. 



IV. SUMMARY AND CONCLUSIONS 

Recently, the CMS and ATLAS collaborations announced the discovery of a new reso- 
nance whose property is in rough agreement with the SM Higgs boson. But the nature of 
this new state, especially its role in electroweak symmetry breaking, needs to be scrutinized. 
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So the most urgent task for the LHC is to test the property of this Higgs-like boson by 
measuring all the possible production and decay channels with high luminosity. Among the 
production channels, the Higgs pair production is a rare process at the LHC. Since it can 
play an important role for testing the Higgs self-couplings, it will be measured at the LHC 
with high luminosity. 

In this work we studied the pair production of the SM-like Higgs boson in the popular 
SUSY models: the MSSM and NMSSM. To make our study realistic, we first scanned the 
parameter space of each model by considering various experimental constraints. Then we 
examined the Higgs pair production in the allowed parameter space. We found that for 
most cases in both models, the dominant contribution to the pair production comes from 
the gluon fusion process with its rate maximized at a moderately light t\ and a large trilinear 
soft breaking parameter A t . The production rate can be sizably enhanced relative to the 
SM prediction: (Tsusy/ctsm can reach 10, and even for a TeV-scale stop it can also exceed 
1.3. For each model we also calculated its x 2 with current Higgs data and found that in the 
most favored parameter region the value of (Jsusy/vsm is approximately 1.45 in the MSSM, 
while in the NMSSM it varies from 0.7 to 2.4. 
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